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Abstract 

We obtain an exact result for the non-perturbative quark (antiquark) production rate and its px 
distribution from a constant SU(3) chromo-electric field E a with arbitary color index a by directly 
evaluating the path integral. Unlike the WKB tunneling result, which depends only on one gauge 
invariant quantity \E\, the strength of the chromo-electric field, we find that the exact result for 
the pt distribution for quark (antiquark) production rate depends on two independent Casimir 
(gauge) invariants, E a E a and [d abc E a E b E c ] 2 . 
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Non-perturbative quark anti-quark pair production from a constant chromo-electric field 
is widely employed to study hadronization at low in high energy e + e~ and pp collisions 
Q]. In these approaches the color flux-tube energy density or the string tension is related 
to the constant chromo-electric field strength \E\. In high energy heavy-ion collisions at the 
RHIC and the LHC l2j a classical chromo field might be formed just after two nuclei pass 
through each other U^J. In order to study the production of a quark-gluon plasma from 
a classical chromo field it is necessary to know how quarks and gluons are formed from the 
latter. In a recent paper \^ we have derived a formula for the rate for non-perturbative 
gluon pair production and its pt distribution from a constant SU(3) chromo-electric field 
with arbitary color index a via vacuum polarization. In this paper we will extend our study 
to quark-antiquark pair production. 

We will not employ Schwinger's proper time method jfj in our calculation. Although 
this method is widely used to obtain the total fermion production rate dN/d 4 x jf| ?| (for 
a review see Jsj]), this method cannot be used to obtain the vt distribution of the rate 
dN / d A xd 2 pT- For this purpose the WKB tunneling method has been used in the past 

to approximate the pt distribution for the quark (antiquark) production rate. Although the 
WKB tunneling result for the pt distribution is correct in QED it is not necessarily true 
in QCD because of the presence of the non-trivial color generators in the fundamental and 
adjoint representations of SU(3). For this reason we will directly evaluate the path integral 
in this paper and obtain an exact result for the pt distribution of the non-perturbative 
quark (antiquark) production rate from a constant chromo-electric field with arbitary color 
in the gauge group SU(3). We find that unlike the WKB tunneling result, which depends 
on one gauge invariant quantity \E\, the strength of the chromo-electric field, 3, [l0| the 
exact result for the pr distribution for the quark (antiquark) production rate depends on 
two independent gauge invariants, E a E a and [d a b c E a E b E c ] 2 . 

We obtain the following formula for the number of non-perturbative quarks (antiquarks) 
produced per unit time, per unit volume and per unit transverse momentum from a given 
constant chromo-electric field E a 

JAT _ 1 3 7r(p2 +m 2 ) 

" "" ~ TTT E A In 1 - e W ], (1) 



dtd 3 xd 2 px 47T 3 , | 

where m is the mass of the quark. This result is gauge invariant because it depends on the 
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following gauge invariant eigenvalues 



Ai = J y cos9 ' 

A 2 = ^ cos (2tt/3 - 6) , 

As = cos (2tt/3 + 0) , (2) 

where 9 is given by 

cos 2 3^ =3C 2 /C 1 3 . (3) 

These eigenvalues only depend on two independent Casimir invariants for SU(3) 

d = £ a £<\ C 2 = [d abc E a E b Ef , (4) 

where a, 6, c = 1,...,8. Note that < cos 2 3# < 1 because Cf > 3C 2 and both C\ and 
C*2 are positive. The integration over in eq. (0) reproduces Schwinger's result for total 
production rate dN/d 4 x 2j- 

can be contrasted with the following formula obtained by the 



The exact result in eq 
WKB tunneling method 



if 



dN - —\oE\ t(?4+™ 2 ) 

sfe = -w- 11111 - e ~^ ] ■ (5) 

In our result in eq. (0) the symmetric tensor d a b c appears. Hence the WKB tunneling 
method does not reproduce the correct result for the pt distribution of the quark (antiquark) 
production rate from a constant chromo-electric field E a . We now present a derivation of 
eq. ©. 

The Lagrangian density for a quark in a non-abelian background field is given by 

£ = ft [(% j> - gT«4 a ) - mSij] ft = ft M i3 [A] ft , (6) 

where p^ is the momentum operator and T?- is the generator in the fundamental representa- 
tion of gauge group SU(3) with a=l,2...8 and i,j = 1,2,3. The vacuum to vacuum transition 
amplitude in the presence of the non-abelian background field A ^ is given by 

j mm f ^ * = = 

/ [#][*] e'/ fc *' M "l°' ** 



The one loop effective action becomes 

S® = - % Tr ln[(^ j> - gT^ a ) - m5 i3 ] + iTi ln[^ j> - mSy] . (8) 

The trace Tr contains an integration over d^x, a sum over color and Lorentz indices and a 
trace over Dirac matrices. Since the trace is invariant under transposition we also get 

S$ = - % Tr ln[(<% $ - gT*4 a ) + m<y + tTr ln[<% j> + mSy] . (9) 

Adding both the above equations we get 

2S$ = - % Tr ln[(<% j> - gT^ f - m 2 b i3 \ + «Tr Info f - m 2 5 lj ] , (10) 

which can be written as 



$ = -iTrH&jP ~ 9T?jA a ) 2 + |^^.F°^-m%]+zTrln[%p 2 -m 2 %]{ll) 
Since it is convenient to work with the trace of the exponential we replace the logarithm 



by 



hi °L = §l[ e i» (b Me) _ e is (a +«)! _ (12) 

b Jos 



We assume that the constant electric field is along the z-axis (the beam direction) and we 
choose the gauge Aq = so that A% = — E a x° [5f. The color indices (a=l,....8) are 
arbitrary. Since Ay = T^E a has three eigenvalues we write after diagonalization 

{k d ) l3 = (Ai, A 2 , A 3 ) . (13) 

The trace over the Dirac matrices (trp) give 

tr D [e is ^ T S F H = 4cosh(^£ a ). (14) 

To reduce this problem to the motion of one harmonic oscillator, we make a similarity 
transformation [a, (we also make a similarity transformation in the group space) and 
obtain 

tr D e isl(5lJ v ~ 9T S Aa ) 2 + § ff M^i^-m 2 <y 
1 cosh(sg(A d ) u )[e ip3po/9( - Ad) e***'-** " m * " 3 2 (^ 2 ) e -i P 3 Po/ 9 (A d ) j y ? ( 15 ) 
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where pr = ypf + p% is the transverse momentum of the quark or antiquark (transverse 
to the electric field direction). Hence from eqs. (fTTj) and (fT2"j) we find 



25$ = i I — E tr 



00 rfs 
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s i=a 



ip 3 Po/g^j g*s(po a -PT 2 - m 2 - 3 2 A 2 aTo 2 +ie) e -ip 3 Po/g\j 



x [4coshs#A.,] - 4e is( P° 2 - p " T2 " m2+ie) ] . (16) 

The trace tr denotes an integral over a complete set of x eigenstates. We add complete sets 
of pj eigenstates, and obtain 



sinhs|gAj| s 

The s-integral at fixed is convergent at s — > 0, but the integration over yields an ex- 
tra factor 1/s so now it seems divergent. However charge renormalization cures this ultravio- 
let problem by subtracting also the term linear in s in the expansion of cosh sgXj/sinh s\g\j\. 
The integral is well behaved as s — > 00. To perform the s contour integration we use the 
well-known expansion 

and then we formally replace s by —is (as first advocated by Schwinger in QED). The 
integral is now real, except for half-circles around the poles at s|gAj| = — inn for n=l,2,3... 
The 1/x term in (jl8j) cancels against the 1/s term in 1)1 7]). This yields the probability for 
quark (antiquark) production per unit time and per unit volume 

W m = 2Im£« = -L /rfV EE" e I5^ ] . (19) 

4?r J n=l i=l n 

Now all that is left is to determine the eigenvalues Aj (j = 1,2,3) of the matrix Ay- = T°-E a 
in the fundamental representation of the gauge group SU(3). Evaluating the traces of Ay, 
A?- and the determinant of A^- we find 

Ai + A 2 + A 3 = 0, (20) 
\l + \ 2 2 + \ 2 3 = l -E a E\ (21) 



and 



AiA 2 A 3 = ^[d abc E a E b E c ], (22) 
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the solution of which is given by eq. (0). 

In this letter we have obtained an exact result for the rate for non-perturbative quark 
(antiquark) production and its p? distribution in a constant chromo-electric field E a with 
arbitary color index a via vacuum polarization. We have used the background field method of 
QCD with the gauge group SU(3). The p T distribution for quark (antiquark) production can 
be applied at the RHIC and the LHC colliders. We find that, unlike the WKB tunneling 
method, the pt distribution of the quark or antiquark production rate depends on two 
independent Casimir (gauge) invariants, E a E a and [d a b c E a E b E c ] 2 . 
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